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Abstract 

> 

OO . We develop a systematic method for solving the 3-loop (^-derivable approx- 



imation to the thermodynamics of the massless cp"^ field theory. The method 
involves expanding sum-integrals in powers of g'^ and m/T, where g is the 



r~| ' coupling constant, m is a variational mass parameter, and T is the temper- 

-L . ature. The problem is reduced to one with the single variational parameter 

0-) : „ 

,^ . m by solving the variational equations order-by-order in g and m/T. At 



the variational point, there are ultraviolet divergences of order g^ that can- 
not be removed by any renormalization of the coupling constant. We define 
a finite thermodynamic potential by truncating at 5**^ order in g and m/T. 
The associated thermodynamic functions seem to be perturbatively stable 
and insensitive to variations in the renormalization scale. 



I. INTRODUCTION 

The thermodynamic functions for massless relativistic field theories at high temperature 
T can be calculated as weak-coupling expansions in the coupling constant g. They have 
been calculated explicitly through order g^ for the massless 0^ field theory [^0, for QED 
0-0, and for nonabelian gauge theories |§-0. Unless the coupling constant is tiny, the 
weak-coupling expansions are poorly convergent and very sensitive to the renormalization 
scale. This makes the weak-coupling expansion essentially useless as a quantitative tool: 
the expansion seems to be reliable only when the coupling constant is so small that the 
corrections to ideal gas behavior are negligibly small. The physical origin of the instability 
seems to be effects associated with screening and quasiparticles. 

A possible solution to this instability problem is to reorganize the weak-coupling ex- 
pansion within a variational framework. A variational approximation can be defined by a 
thermodynamic potential Q that depends on a set of variational parameters rrii. The free 
energy and other thermodynamic functions are given by the values of Q and its derivatives 
at the variational point where dQ/dnii = 0. A variational approximation is systematically 
improvable if there is a sequence of successive approximations to Q that reproduce the weak- 
coupling expansions of the thermodynamic functions to successively higher orders in g. One 
example of a systematically improvable variational approximation is screened perturbation 
theory, which involves a single variational mass parameter ^. 

A variational approach can be useful only if the correct physics can be captured by 
appropriate choices of the variational parameters. Information about screening and quasi- 
particle effects is contained within the exact propagator of the field theory. The possibility 
that these effects are responsible for the instability of the weak-coupling expansion suggests 
the use of the propagator as a variational function. Such a variational formulation was con- 
structed for nonrelativistic fermions by Luttinger and Ward |T^ and by Baym |ll|] around 



1960. It was generalized to relativistic field theories by Cornwall, Jackiw and Tomboulis ||12|| . 
In the case of a relativistic scalar field theory, the propagator has the form [P^ + n(P)]~^, 



where n(P) is the self-energy which depends on the momentum P. The thermodynamic 
potential has the form 

fio[n] = Ig [log (p^ + n) - -p^] + $[n] , (i) 

where the interaction functional $[n] can be expressed as a sum of 2-particle- irreducible 
diagrams. It is constructed so that the solution to the variational equation 6Q0/6II = is 
the exact self-energy, and the value of Qq at the variational point is the exact free energy. 
We can obtain a systematically improvable variational approximation by truncating $ at 
n'th order in the loop expansion, where n = 2, 3, .... We refer to this approximation as the 
ra-loop ^-derivable approximation. 

One remarkable feature of the 2-loop ^-derivable approximation is that the expression 
for the entropy reduces at the variational point to the 1-loop expression. This was shown for 
QED by Vanderheyden and Baym [|13[. It was generalized to QCD by Blaizot, lancu, and 
Rebhan |jT^, who have used it as the basis for a quasiparticle model for the thermodynam- 
ics of the quark-gluon plasma. An alternative quasiparticle model motivated by the 2-loop 
^-derivable approximation for QCD has been developed by Peshier [jT5[. One disadvantage 



of these models is that they are based on special properties of the 2-loop ^-derivable ap- 
proximation. Thus unlike the ^-derivable approximation itself, they are not systematically 
improvable. 

^-derivable approximations are conserving approximations, which means that they are 
consistent with the conservation laws that follow from Noether's theorem. Baym showed 



that an approximation is conserving if and only if it is ^-derivable for some functional $ [11 



The fact that they are conserving may make <l>-derivable approximations particularly useful 



for nonequilibrium problems [16|. They have already proved to be useful for nonequilibrium 



problems in 1 -|- 1 dimensions |T^ 



While the ^-derivable approximation is easily formulated, it is not so easy to solve. If 
the solution for the self-energy 11 (P) is independent of the momentum P, the <l>-derivable 
approximation can be reduced to a single-parameter variational problem that is easily solved 



numerically. An example is the 2-loop ^-derivable approximation for the massless 0'* field 
theory |18[. If the solution for n(P) depends on P, the variational equation is a nontrivial 
integral equation. What makes it complicated to solve is that it is really a coupled set of inte- 
gral equations for the function Il{27mT, |p|) corresponding to all the Matsubara frequencies 
27mT, 72 = 0, ±1, ±2, .... An even more severe obstacle is that the thermodynamic potential 
has ultraviolet divergences that vanish at the variational point only if $ is calculated to all 
orders. They do not vanish away from the variational point, and they do not vanish at the 
variational point if the loop expansion for $ is truncated. 

In this paper, we solve the 3-loop ^-derivable approximation for a massless scalar field 
theory with a (p^ interaction. In section II, we define the coupling constant for the massless 
0^ field theory and give the weak-coupling expansions for several thermodynamic quantities. 
In section III, we summarize the basic features of the ^-derivable approach. In section IV, we 
solve the 2-loop ^-derivable approximation. We show that if the thermodynamic potential 
is expressed in terms of the true coupling constant, it has ultraviolet divergences of order 
g^. In section V, we solve the 3-loop ^-derivable approximation by expanding sum-integrals 
systematically in powers of g and m/T. We show that the thermodynamic potential has 
ultraviolet divergences of order g^. We construct a finite thermodynamic potential by adding 
a term proportional to the square of the variational equation and then truncating after terms 
of 5*^^ order in g and m/T. It defines a stable approximation to the thermodynamic functions 
that is rather insensitive to the choice of renormalization scale. In the concluding section, 
we discuss the outlook for extending our methods to gauge theories. 

II. DEFINING THE THEORY 

We would like to solve the ^-derivable approximation to the thermodynamics of a mass- 
less scalar field theory with a 0^ interaction. The theory can be defined by choosing a 
regularization scheme and specifying the value of the bare coupling constant go. The la- 
grangian for the theory is 
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We choose to use dimensional regularization in 4 — 2e space-time dimensions. The bare 
couphng constant qq then has dimensions (mass)*^. We define the renormahzed couphng 
constant g{^) by the modified minimal subtraction (MS) prescription with renormalization 
scale /i. For convenience, we also define a = (7^/(47r)^ and cto = Qq/IAtt)'^. The relation 
between the bare and renormahzed coupling constants is 
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where ({z) is the Riemann zeta function. In the MS scheme, there is usually also a factor 
of (47re~''')^^ on the left side. We prefer to absorb this factor into the measure of the dimen- 
sionally regularized integrals. The beta function for the running coupling constant a{fi) has 



been calculated to five- loop order |T9|. The first few terms can be obtained by differentiating 
(I) with respect to /i and using the fact that ao is independent of fi: 



fj.—-a = —2ea + P{a) 



(4) 



where the beta function is 



/3(a) = 3«2 _ 11^3 ^ j^ll^ ^ i2C(3) ) a^ + 0{a'' 
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The thermodynamic functions for this theory are known to order g^ [0,|l],3- The weak- 
coupling expansion for the free energy density is 



J^ = ^i 



+540 1^2 log ^ - 3L - ^ + 4 log 2 - 7 + '^^^] (f ) '^' + 0{a' log a) 



(6) 



where a = a{fi), L = log(/i/47rT), 7 is Euler's gamma constant, and J^ideai is the pressure 
of the ideal gas of a free massless boson: 
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The weak-coupling expansions for the pressure, entropy density, and energy density are given 
by: 
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where the partial derivatives in (^) and (|10D are evaluated with a(/i) and /i held fixed. The 
free energy density and other thermodynamic functions satisfy simple renormalization group 
equations: 






J^ = 0. 



(11) 



If the power series expansion for (3 (a) and the weak-coupling expansion for JF are truncated 
at some order in a, then the renormalization group equation (^) is satisfied only up to 
higher orders in a. 

The renormalization group can be used to improve the weak-coupling expansions of 
the thermodynamic functions. This can be accomplished for the free energy (^) simply by 
setting /i = c{2ttT), where c is a constant, in the logarithm L and in a(/i). In this case, 
there are alternative thermodynamic definitions of the entropy and energy density in which 
the partial derivatives in @ and (|l^) are replaced by total derivatives with respect to T. 
For example, the thermodynamic entropy is defined by 



d 
It differs from the entropy (j^) by a term proportional to the left side of ([Tl|): 



(12) 



S, 



th 



s = _1^ 
fidT 



9 n, \ d 



T. 



(13) 



The right side is nonzero if the expansions of (3 and T in powers of a are truncated. 

The poor convergence of the weak-coupling expansion is illustrated in Figure 1, which 
shows the free energy divided by that of the ideal gas as a function of g{2TiT). The solid 
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FIG. 1. Weak-coupling expansion for the pressure divided by that of the ideal gas as a function 
of g{2TTT). Solid lines correspond to the truncation at order g"', n = 2, 3, 4, 5 for /i = 2ttT. Dashed 
lines correspond to the truncation at order g^ for ^ = ttT and fi = AttT. 

curves are the weak-coupling expansion with renormalization scale fiT = SttT truncated 
after orders g"', n = 2,3,4,5. They show no sign of converging except at extremely small 
values of g. The sensitivity of the weak-coupling expansion to the renormalization scale 
is also illustrated in Figure 1. The dashed lines are the weak-coupling expansion through 
order g^{fi) for fi = irT, fi = AnT plotted as a function of gifJ^r)- The relation between g{fi) 
and g{^T) is obtained by integrating numerically from 27rT to /i using the beta function (^ 
truncated after the a^ term. Note that g = A, which corresponds to a = 0.1013, represents 
relatively strong coupling in the sense that it is close to the Landau pole. If g{27cT) = 4, 
the Landau pole at which the running coupling constant defined by the 3-loop beta function 
diverges is /i = 12.15(27rT). 

Another observable that has been calculated to relatively high order is the screening 
mass p|. It is known up to corrections of order g^: 
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FIG. 2. Weak-coupling expansion for the screening mass divided by the leading order result 
mf^Q = (7^(27rT)r^/24 as a function of g{2-KT). Solid lines correspond to the truncation at order 
g^, n = 2, 3, 4 for ^ = 2'kT. Dashed lines correspond to the truncation at order g'^ for /i = vrT and 
// = 47rr. 

The convergence of the weak-coupling expansion is illustrated in Figure 2, which shows 
m^ divided by the leading order expression mLo = 5'^(a^t)7'^/24 as a function of g{27iT). 
The solid lines are the weak-coupling expansion ( [T^ ) truncated after the term of order g"', 
n = 2, 3, 4. The weak-coupling expansions to orders g^ and g"^ agree only at extremely small 
values of g. The expansion to order g'^ is also very sensitive to the renormalization scale as 
illustrated in Figure 2. The dashed lines show ml/m^Q, with ml truncated after the g'^^n) 
term as a function of g{2TTT) for /i = vrT and /i = AttT. 

III. ^-DERIVABLE APPROXIMATIONS 



^-derivable approximations are based on the skeleton expansion of the free energy [10 



which expresses it in terms of exact propagators. The exact propagator D{P) for a massless 
scalar theory has the form 



The skeleton expansion for the free energy defines a functional Q[D] of the propagator 
that we will refer to as the thermodynamic potential. The exact propagator satisfies the 
variational equation 6Q/6D = 0. The free energy is therefore the variational minimum of 
the functional r2[D]: 



J^ = n[D] 

sn/5D=o 
The thermodynamic potential can be expressed in the form 



(16) 



n[D] = \g log (p^ + n(p)) - Ig n(p)-^,^^L_ + ^n], (17) 

where the functional $[-D] contains all the effects of interactions. It can be expressed as a 
sum of 2-particle-irreducible Feynman diagrams. The variational equation 6Q/6D = can 
be written 

The variational derivative is calculated with the bare coupling constant Oq and the ultraviolet 
cutoff parameter e held fixed or, equivalently, with the renormalized coupling constant a and 
the renormalization scale /z held fixed. 

An approximation is called ^-derivable if it can be obtained by making the same approx- 
imation for the functional $ in both (|l^ and (|1^) [|ll|]. ^-derivable approximations have 
several attractive features. One such feature is self- consistency: the thermodynamic poten- 
tial ([T7|), which includes the effects of interactions through the $ term, is evaluated with a 
propagator that takes into account those same interactions through the variational equation 
(|18|). Another attractive feature is thermodynamic consistency. The entropy density can be 
expressed as 
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where d/dT refers to the partial derivative with D{P), ao, and e held fixed. The thermo- 
dynamic entropy can be written 



(20) 

5n/5D=0/ 



where the derivative is evaluated with ao and e fixed. The equivalence of (|T^ and 
follows from the variational equation 6Q/6D = 0. 

The n-loop ^-derivable approximations defined by truncating $ at successively higher 
orders in the bare coupling constant, or equivalently the loop expansion, are particularly 
useful. They define a systematically improvable variational approximation. If $ is truncated 
at order g^ , the free energy will after renormalization include correctly all terms in the 
weak-coupling expansion through order f?^""^, as well as an infinite series of higher terms 
in g. Thus each successive approximation will have an error that is suppressed paramet- 
rically by one higher order in g"^. The propagator for each approximation is given by the 
self-consistent solution to a variational equation that takes into account the screening and 
quasiparticle effects characteristic of a plasma. Thus if the poor convergence of the weak- 
coupling expansion is due to treating plasma effects as a perturbation, the n-loop ^-derivable 
approximations could provide the solution to this problem. 

Unfortunately, ^-derivable approximations have a problem that has limited their use- 
fulness in relativistic field theories. The problem is that the thermodynamic potential fl 
is generally ultraviolet divergent. There are two classes of ultraviolet divergences. One 
class consists of ultraviolet divergences that vanish at the variational point. They provide 
a serious complication for devising numerical methods for minimizing the thermodynamical 
potential. The other class of ultraviolet divergences are those that survive at the variational 
point. They arise from having truncated $ at n*^ order in the loop expansion. Renor- 
malization theory guarantees that the weak-coupling expansion for the free energy will be 
finite through order g^^^^, but there will generally be ultraviolet divergences at order g^^. 
In the 2-loop ^-derivable approximation for the massless 0^ field theory, these ultraviolet 
divergences can be hidden by introducing a coupling constant that has the wrong beta func- 
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tion. However in the 3-loop ^-derivable approximation, these ultraviolet divergences must 
be dealt with. 

IV. TWO-LOOP ^-DERIVABLE APPROXIMATION 

In the 2-loop ^-derivable approximation, the solution for the exact self-energy is inde- 
pendent of the momentum. We anticipate this fact by denoting the exact self-energy by 
m'^. The variable m can be interpreted as the screening mass. The solution to the 2-loop 
<l>-derivable approximation has been given previously in Refs. |l^ and |l^. We first solve 



the 2-loop ^-derivable approximation exactly in terms of a coupling constant g{^) that has 
the wrong beta function. We show that the expansion of the thermodynamic potential in 
powers of m/T converges much more rapidly to the exact answer than the weak-coupling 
expansion in powers of g. We then show that when the solution is expressed in terms of the 
true coupling constant g, there are unavoidable ultraviolet divergences at order g^. 
The thermodynamic potential in the 2-loop ^-derivable approximation is 

A''Oo(m) = ipo^iP^ + ."') - I'n^i J^ + l(9o,-r {g ^)' , (21) 

where go is the bare coupling constant. Our definition of the dimensionally regularized sum- 
integral is given in ( p.l|) and includes a factor of fi'^'^ in the measure. Each term in (pTf) has 



an overall multiplicative factor of /i^'^, and therefore fio(m) is independent of /i for fixed go 
and e. This factor of fi'^'^ is convenient, because it gives every term in (PH) the dimensions 
(mass)"^, even for e 7^ 0. The variational equation dQ/dm = reduces to the simple gap 
equation 



m^ 
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The 2-loop ^-derivable approximation for the free energy density JF is obtained by inserting 
the solution to the gap equation (|2^) into the thermodynamic potential (PT|). 
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A. Finite thermodynamic potential 



The sum-integrals that appear in the thermodynamic potential ( pTj) and in the gap 
equation (^) are ultraviolet divergent. Ideally, we would like to have a thermodynamic 
potential i7(r7i) that is a finite function of the variational parameter m. To construct such 
a function, we begin by making the ultraviolet divergences explicit as poles in e. The sum- 
integrals can be expressed in the form 
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where m = m/{2nT), ft = /i/(27rT), and a(e) and 6(e) are constants that approach 1 in the 
limit e — i> 0: 
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The function /(m^, e) is given by the integral 
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and f'{rh%e) is its derivative with respect to rh . In the limit m — * 0, f{rh,0) approaches 
m^. The sum-integral (E^ can be obtained simply by differentiating (ESf) with respect to 



m 



^. Inserting ( PBj ) and (^) into (pij), the thermodynamic potential reduces to 
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The gap equation (|^) reduces to 
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We wish to construct a finite thermodynamic potential Q that has the same value at the 
variational point as (|28|). The gap equation (|29| ) can be expressed in the form rh^—G = 0. We 
can add any multiple of (m^ — G)^ to ([28| ) without changing its value at the variational point. 
We can use this freedom to cancel the double pole in e in (p8|). Adding 45(m^ — G')^/(Q;o/i~^'') 
to (PBD, the resulting function reduces to 
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This thermodynamic potential will be finite as e — i> provided that the "renormalized 

coupling constant" a{fj.) defined by 

1 11 

(31) 



is finite as e ^ 0. This equation does not define the value of a(yu) at any specific scale fi, 



because (|3T1) does not have a smooth limit as e — > 0. However its dependence on /x is well 



defined. Differentiating (|3T|) with respect to /i, we obtain 



/i-— a(/i) = — 2ea(/i) + a^(/i) . 



(32) 



This equation has a smooth limit as e — > 0, so the /i-dependence of a{fi) is well-defined. Note 
that the coefficient of a^ in (^) differs from the coefficient of a^ in the true beta function 
(^) by a factor of 3. 

Using (^H) to eliminate ao from ( ^0|) in favor of Q;(/i), we obtain a finite thermodynamic 
potential in the limit e — *> 0: 
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Upon varying (|33D with respect to m^, we obtain the gap equation 
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The free energy !F(T) in the 2-loop ^-derivable approximation is obtained by solving (|35D 
for m and inserting the solution into (p^). 

The thermodynamic potential ( P^D is independent of the scale /i, because it satisfies the 
renormalization group equation 



djji da 



n{m) = . (36) 



However the ^-derivable free energy still has an ambiguity associated with the renormaliza- 
tion scale. The ambiguity arises because ( pi]) determines the dependence of Q;(yu) on // but 
it does not define the value of a(/i) at any scale /i. The value of a(/i) becomes well-defined 
only after some initial value a{no) is specified. A convenient choice for /iq is the scale at 
which a has the same value as the true coupling constant a: 

a{fXo) = a{fXo) ■ (37) 

Since a and a have different beta functions, their values are different at other values of /i: 

a^^) + 0{a^). (38) 



a{n) = Q;(/i) — 2 log — a^(/i) + 
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The ambiguity in the matching scale /xq is the price that must be paid for using a coupling 
constant with the wrong beta function. A reasonable choice for the matching scale is /xq = 
27iT. However in estimating the errors of the ^-derivable approximation, one should take 
into account the ambiguity in /io, perhaps by allowing it to vary by a factor of 2 around the 
preferred value. 

B. Weak-coupling Expansion 

The two-loop free energy in the <l>-derivable approximation can be expanded systemati- 
cally as a weak-coupling expansion in powers of a^'^. One can obtain a sequence of analytic 
approximation to the free energy by truncating the expansion. Unfortunately we will find 
that the series converges rather slowly. 
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To express the free energy as an expansion in powers of a^/^, we need to expand the 
solution to the gap equation (|35|) in powers of a^/^ and insert it into (|33|) . We first need to 



expand the function f{rh?) in powers of m. The expansion, which is derived in Appendix 0, 
has the form 

f{rh^) =rh^^Am^-?, (log - - - + 7 j m^ + ^ f^rh^- , (39) 

where the coefficients /„ are 

^ 6r(n-|)C(2n-3) 

•'" * ' r(i)n! ■ ^ ' 

Note that the only odd power of rh is the m^ term. The asymptotic behavior of the coeffi- 
cients as n — > cxD is /„ ^ (— l)"'+-'^(6/v^)'^~^^^- Using the ratio test, we can infer that the 
region of convergence of the power series in P5| ) is |?Ti| < 2txT. 

Inserting the expansion (|39D into (^), the thermodynamic potential becomes 



^^""^ -l_15|m2-4m3-3f- + L + 7)m^ + ^/„m2"l , (41) 
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where L = log(/i/47rT). Differentiating the expansion (|39|) and inserting it into (^), the 
gap equation becomes 

m^ = - 1 1 - 6m - 6(L + 7)771^ + X] nf^m^''-^ I . (42) 

" I n=3 J 

Using the ratio test, we can infer that the regions of convergence of the series in (^) and 



(^) is \m\ < 2ttT. As long as m is inside this radius of convergence, the functions of m/T on 
the right can be evaluated simply by adding terms in the series until the desired numerical 
accuracy is achieved. The first few terms in the series for (|43) are 



2 « 



m- = - {1 - 6m - 6(L + 7)771^ + 1.803m^ - 0.778m^ + ...}. (43) 

For asymptotically weak coupling, the solution is m = gT j\f2A. The radius of convergence 
in m translates into |^| < 47r-\/6 '^^ 30.8. 

The solution to the gap equation (^21) up to errors of order ci^l'^ is 
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(a) (b) 

FIG. 3. 2-loop <I>-derivable approximation to the pressure divided by that of the ideal gas (heavy 
sohd hne) as a function of g{2TrT). bi Figure (a), this is compared to the weak-couphng expansion 
truncated at orders g", n = 2,3,4,5,6. In Figure (b), it is compared to the m/T expansion 
truncated at orders 5", n = 2, 3, 4, 6. 
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where a = a{fi) and L = log(/i/47rT). The weak-coupling expansion for the two-loop $- 
derivable free energy is obtained by inserting the expansion (|4^) into (^Tj) and expanding in 
powers of a^/^. The expansion up to errors of order a"^/^ is 
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The convergence behavior of the weak-coupling expansion for the 2-loop ^-derivable 
approximation is illustrated in Figure 3, which shows the free energy divided by that of the 
ideal gas as a function of g{2TxT). The exact result is shown as a heavy solid line. The 
light solid lines show the weak-coupling expansions for /i = 2tiT truncated after the terms 
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of order g"', n = 2,3,4,5. The weak-coupling expansion converges to the exact answer, 
which is reflected in the fact that successive approximations track the exact answer out to 
higher values of g. However the convergence is relatively slow. While the exact answer is 
completely independent of the renormalization scale /i, the weak-coupling expansion has a 
strong sensitivity to /i as illustrated in Figure 3a. The dashed lines are the weak-coupling 
expansions for /i = vrT and /i = AnT truncated after the g^{fi) term as a function of g(2TxT). 
The relation between g{^) and g{2TiT) is given by the exact solution to the renormalization 
group equation: 

aiii) = z ^ '—-, . (46) 

^^' 1 - a(27rT) log(/i/27rT) ^ ^ 

Note that the Landau pole appears at /i = e^^°'^'^^'^\27rT). For g{27rT) = 8, it appears at 
/i = 11.79(27rT). Thus g = 8 in the 2-loop ^-derivable approximation is roughly as close to 
the Landau pole as g = 4, where g is the running coupling constant defined by the 3-loop 
beta function of the 0'^ theory. 

C. Expansion in m/T 

The slow convergence of the weak-coupling expansion provides motivation for developing 
a better sequence of successive approximations to the ^-derivable free energy. We would like 
an approximation scheme that is applicable not only to the two-loop ^-derivable approx- 
imation, but can be extended to higher loops. We will develop an approximation scheme 
that is based on expanding sum-integrals in powers of m/T. 

The expansion of the finite thermodynamic potential Q{m) in powers of m = m/{2TiT) 
is given in (0). We can define a series of successive approximations Vt^"'\m) by truncating 
after terms of n-th order in f^, with m treated as order g. For example, the thermodynamic 
potential truncated after the term of order g^ is 

— ^ ' - 1 = — - m^ + 4m=^ + 3 (L + 7) m^ - -C,{^)m^ . (47) 

15 y T\Ac&\ jo. 2 
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The corresponding gap equation is 

(48) 



- 2 (^ 

m = — 
6 



3 

1 — 6m — 6(L + 7)m^ H — C(3)m^ 



The first few of these successive approximations to the free energy can be calculated 
analytically. The n = 2 truncation includes only the first two terms of the right side of (^^ . 



The solution to the gap equation is simply m = J a/6. The resulting expression for the free 
energy is 

^ = l--a. (49) 

-'ideal ^ 

The n = 3 truncation includes also the rh^ term in (|^). The gap equation is a quadratic 
equation whose solution is 



,6 4 2 

The resulting expression for the free energy is 



1 _ 1_ 1 _ 
m = \j-a + -a'^ a. (50) 



^(3) 5 
1 a 



a a^ 



l + 6a + 6a^-4(2 + 3a)W-+ ^ 



(51) 



•'ideal ^ 

In the strong coupling limit a — i> oo, the solutions for the n = 3 truncation have finite 
limits: m — > | and J^'-^V-^ideai -^ %■ Note that solutions ( ^Ol) and dSTf ) both depend on the 



renormalization scale /i through a{^). The screening mass and the free energy JF^^) for the 
n = 4 truncation can be obtained simply by substituting a — > 1/(1/q; + L + 7) into (^) and 
(pT|). The resulting expressions are completely independent of /i. 

The thermodynamic functions defined by truncating the m/T expansion converge much 
more quickly than the weak-coupling expansion. This is illustrated in Fig. 3b, which shows 
the free energy divided by the ideal gas as a function of g{2TtT). The heavy solid line is the 
exact 2-loop ^-derivable result. The slighter solid lines are the results from truncating the 
thermodynamic potential after terms of order (yf", n = 2,3,4,6, with m treated as order g. 
For n = 2, 3, we set /i = 27rT, but beginning at order g'^, the result is completely independent 
of yU. The results for n = 4, 6 are so close to the exact result that they cannot be distinguished 
from it in the Figure. Thus the convergence of the m/T expansion is remarkably fast. 



D. Solution in terms of True Coupling Constant 

In the 3-loop ^-derivable approximation, we will not have the luxury of being able to 
calculate the thermodynamic potential as explicitly as in (pHI). The best we will be able to 
do is calculate the sum- integrals as truncated expansions in m/T. We will also be unable 
to absorb all the ultraviolet divergences in the free energy into a renormalized coupling 
constant a{^) that runs differently from the true coupling constant a{^). We will therefore 
repeat the solution to the 2-loop ^-derivable approximation in a manner that is as parallel 
as possible to the method we will use for the 3-loop case. 

Our goal will be to construct a thermodynamic potential that depends on the true 
coupling constant a{n) defined by (|) and is finite to as high an order in g as is possible. 
We have the freedom to remove ultraviolet divergences from the thermodynamic potential 
(pl|) without changing the free energy by adding multiples of {m? — G)^, where G is the 



right side of the gap equation (p2|). We can cancel the divergences proportional to f^g/e^ in 
the last term of (|2T|) by subtracting {m? — GY / {2{gQii~''Y) . The resulting thermodynamic 
potential is 

^^^^n{m) = \i log(p2 + m") - "^1 . (52) 

2Tp 2{golJ, ^Y 

Expanding in powers of m/T, we get 

^ll^}^ = 1 _ 15 1^2 _ 4^3 _ 3 l^^^ + 1 + L + 7^ m4 + 0{m^) . (53) 

^idcai i V"oyU ^' 2e ) ] 

If we substitute the expression (|^) for the bare coupling constant in terms of the true 
renormalized coupling constant a(/i), we find that there are still ultraviolet divergences 
proportional to m^: 

^^!^}!^ = 1 _ 15 (^2 _ 4^3 _ 3 [1 + ('.i + 2. + 7) - —a\ m4 + 0{g^) . (54) 
J^ideai l La V e / 12e J J 

The closest we can come to constructing a finite thermodynamic potential that is a function 
of the true coupling constant is to truncate after the terms of order g^: 
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The corresponding gap equation is 

m^ = — [1 — 6m] . (56) 

Note that these expressions differ from those obtained by truncating (^) and (^), because 
they involve the true couphng constant a instead a. The solution to the gap equation ( |5^ 
is obtained by substituting a for a in (|50|) . The resulting expression for the free energy is 
obtained by substituting a for a in (0). 

V. THREE-LOOP ^-DERIVABLE APPROXIMATION 

In the three-loop ^-derivable approximation, the solution to the variational equation (piSD 
for the self-energy 11 (P) is a nontrivial function of the momentum P. The thermodynamic 
potential in the 3-loop ^-derivable approximation is 

where 53baii is the basketball sum-integral: 

^•""^ " TpQR p2 + n(P) Q2 + n(Q) i?2 + n(i?) 52 + n(5) ' ^^^^ 

where S = —{P + Q + R). The variational equation obtained by varying (|57D with respect 
to n(P) is 

n(P) = l(,.,-')% ^.^ - i(*.-T3...(P) . (59) 

where Q'sun(P) is the sunset sum-integral: 

Ssun(P) = f^^ Q2 + n(Q) p2 + n(p) 52 + n(5) ' ^^°^ 

The three-loop ^-derivable approximation to the free energy is obtained by solving ( |5D| ) for 
n(P) and inserting the solution into the thermodynamic potential 
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A. Mass Variable 

Our strategy for solving the 3-loop ^-derivable approximation is to introduce a mass 
variable m which is of order gT in the weak-coupling limit and then calculate the sum- 
integrals in (|57| ) and ([59| ) as double expansions in ao and m/T. We will find that the gap 
equation (|^) for n(P) has a recursive structure that allows us to solve for its dependence 
on P. Inserting the solution into (|57D, we can reduce the thermodynamic potential Q{m) to 
a double expansion in oq and m/T. 

The simplest choice for the variational parameter conceptually is the self-energy at zero 
external momentum 11(0). This is simply the value of the variational function n(P) at a 
particular point in momentum space. One of the disadvantages of n(0) as a variational 
parameter is that it is ultraviolet divergent. Renormalization theory guarantees that in 
conventional perturbation theory the ultraviolet divergences in the exact propagator 1/(P^-|- 
n(P)) are canceled order-by-order in the renormalized coupling constant a by multiplying 
by Z^^, where Z is the wavefunction renormalization constant. Thus Zll{0) should be finite 
order-by-order in a. Since Z is divergent at order a^ and n(0) begins at order a, n(0) 
must have ultraviolet divergences at order a^. Now the ^-derivable approximation involves 
a truncation of perturbation theory, so multiplicative renormalizability is guaranteed only 
up to the truncation error. Since the three-loop ^-derivable approximation does not include 
all three-loop self-energy diagrams, the truncation error for Zn(0) is of order a^. If we use 
n(0) as a variational parameter, ultraviolet divergences in the solution to the gap equation 
will come partly from true divergences of 11(0) and partly from the failure of multiplicative 
renormalizability in the 3-loop ^-derivable approximation. 

A choice for the variational parameter with more convenient ultraviolet properties is the 
screening mass defined by 

n(0, p) — *> m^ a.sp'^^—ml. (61) 

Renormalization theory guarantees that in the conventional perturbation theory ml is finite 
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order-by-order in the renormalized coupling constant a. Since the three-loop ^-derivable 
approximation does not include all three-loop self-energy diagrams, m^ is guaranteed to be 
finite only up to corrections of order a^ . But any divergences at this order will come from 
the failure of multiplicative renormalizability in the 3-loop ^-derivable approximation. 

In the three-loop ^-derivable approximation, the screening mass rus defined by (^Tj) is 
the solution to the gap equation 

m^ = \^9.^-')% Q. /n(g) - ^(^0/."^)'^™ , (62) 

where we have introduced the short-hand 



^screen — ^sunl^U, PJ 



(63) 



p=im 



Since the sum- integrals in (|62D can be expressed as expansions in 11(0) /T, this equation 
defines rn^ as a function of n(0). We will find it more convenient to calculate the sum- 
integrals as expansions in m/T. Subtracting (|6^ from (|59D , the variational equation reduces 
to 

n(P) = m^ - ligofl-')' [%un{P) - ^screen] • (64) 

If we set P = and evaluate the sunset sum-integrals on the right-side as expansions in 
m/T, we get an equation for 11(0) as a function oim. This guarantees that the mass variable 
m defined by (|6l|) can be used as a variational parameter in place of 11(0). The solution to 
the gap equation (|62D for m will be the ^-derivable approximation to the screening mass 
rris. 

B. Solution to the Variational Equation 



In the variational equation ( |6^ for the self-energy n(P), the sum-integrals extend over 
all momenta. There are two important momentum scales: the hard scale 2ttT and the soft 
scale m. The hard region for the momentum P = {27mT, p) includes n ^ for all p and 
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also n = with p of order T. The soft region is n = and p of order m. We will solve the 
variational equation in the two momentum regions separately. 

In order to solve the variational equation for n(P), we will assume that the scales m and 
27rT can be separated. This will allow n(P) to be expanded in powers of g^^^^ and m/T. 
In the variational equation (|5^, the bare coupling constant appears only in the combination 
((7oAi~^)^- We therefore assume that the self-energy n(P) can be expanded in powers of 
{gQ^~^)'^, with coefficients that are functions of P, m and T. These functions in turn can be 
expanded in powers of m/T. For hard momentum P, we write the expansion in the form 

n(P) = m' + {gofi-T [n4,o(P) + n4,i(P) + n4,2(P) + n4,3(P) + . . .] 

+ (^o/i-^)' [n8,-2(P) + n8,_i(P) + ...] + ... for hard P, (65) 

where n„ ^(P) is of order T^(r?7,/T)^ when P is of order T. We have anticipated the solution 
to the variational equation by writing explicitly all terms that contribute through order 
g'^T'^. For soft momentum P = (0, p), we write the self-energy as 11(0, p) = m'^ + cr{p). The 
m/T expansion of the function a{p) has the form 

o-(p) = (fi'oyu"')*^ [o"4 -2(p) + o-4,o(p) + o-4,i(p) + • • •] 

+ (goi^-'f [asMp) + ■■■] + ••• for soft p , (66) 

where (Jn,k{p) is of order m'^{m/TY when p is of order m. We have anticipated the solution 
to the variational equation by writing explicitly all terms that contribute through order 
g^m?. Upon inserting the expansions (|65|) and (^) into the variational equation (0) and 



expanding systematically in powers of gofJ'~'^ and m/T, we will find that it has a recursive 
structure that either determines n„^j!c(P) and (Jn,k{p) or allows them to be expressed in terms 
of lower order functions. 

The m,/T expansions for the sunset sum- integral Q'sun(-P) are derived in Appendix ^. 
For hard momentum P, the m/T expansion including all terms up to errors of order g^T'^ 
is the sum of ( lAl^ ), ( lAlil ), and (lAlel) : 
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P2 rQR{Q^yR^{P + Q + R) 



The momentum integral Ji is proportional to m^"^'^ and is given in ( |C.7| ). We will find 
that the momentum integral involving cr4,_2 is proportional to m~^~'^'^. For soft momentum 
P = (0,p), the sunset sum-integral, including all terms up to errors of order g'^T'^, is the 
sum of { \K7[Ei ), ( ^:^ ), and 1^;^): 



53su„(0, p) = T^Isnnip) + 3Th^ -^ 



QR (Q2)2i?2(Q + ^)2 ■ ^^ JPq^ (g2)3i?2(Q + ^)2 

2 

e\4 



. a4,-2(r) /: 1 2 /• ^4,-2 (g) Hn + nh 

/ V2 + ^2)2J(^(g2)2+^ 7q(^2 + ^2)2^bubUP + qU 



+0{g^T^) . (68) 



The sunset integral Isnn{p) is defined in ( |C.10|) and scales like ?ti ^ ^^ when p is of order m. 
The bubble integral Ihuh{p) is defined in ( |C.9| ) and scales like m~^^^^ when p is of order m. 
We proceed to solve the variational equation ( |64D for n(P). We first consider the equa- 
tion for hard momentum P. On the left side of (0), we insert the expansion ( ]65| ) for 
n(P). On the right side of (p^, we insert the m/T expansions (p7|) and (p5|) of the sunset 



sum-integrals. Matching coefficients of ((^oA* ^)^" on both sides of the equation and then 
identifying the terms in the m/T expansions, we find 



2 ^P2 ' 2 J^iJ (Q2)2/?2(p + g + ^)2 
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{2/dW 



n« 



o^4-2(g) 



\rp2 

2 Jdiq'^ + my 



/bub(|k + q|) 



k=im 



(71) 
(72) 



We next consider the variational equation for soft momentum P = (0,p). On the left side 
of (0), we insert the expansion (BBp for 11(0, p). On the right side of (p3), we insert the 



m/T expansions (|68D of the sunset sum-integral at soft momentum. Matching coefficients of 
{qqU'^Y'^ on both sides of the equation and then identifying the terms in the m/T expansions, 
we find 

1, 



D 



{A/dW 



a8,_4(p) = -T'^ / 71^T^-^\2^bub(|P + q|) + n8,_2 . 

2 Jq (g^ + ni'^Y 



(73) 
(74) 

(75) 



Notice that the matching equations in (§9|)-(^) have a recursive structure. The functions 
1140, n4 1, 114 2, o"4,_2, and a^^ are given explicitly in terms of integrals and sum-integrals. 
The functions Ilg _2 and 0"8^_4 are given in terms of integrals of the function 0"4^_2(q') which 
has already been determined. The recursive structure of these equations continues at higher 
order in g. Thus, by using this method, we can solve the variational equation for n(P) 
order-by-order in g^ and m/T for both hard P and soft P. 

In order to calculate the thermodynamic potential, we will need to calculate the sum- 
integrals and integrals involving some of the functions Yim,k{,P) and am,k{'p)- Using the 
solutions (^) and ([70D, we obtain the following sum-integrals: 

n4,o(P) U 1 1. 



2^2 



p {P 

n4,i(P) 
p {P 



QfpQR{P^yQ^R^{P + Q + Rf 6 



T^'I 



2\2 
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{im^ 
2' 



p{P 



2^2 



fpQ {P^YQ^P + QY \TP (P2)2 
Using the solutions ([731) -([75|), we obtain the following integrals: 

0-4 -2(p) I r^3 f Isnnip) - Isnn{im) 



(76) 
(77) 
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(78) 



a4,ob) ^rr^rJ^ f 1 (4/^)P' 



6 'fpQ WP2)2Q2^P , Qp ^p2^3Q2^p , Qp J ' ^ ) 



p(p2 + m2)2 6 5^Q V(^')'Q'(^ + Q)' (P2)3g2(p + g)5 



Vp (p2 + m2)3 ~ 36 Jp (p2 + m2)3 ' ^ ^ 

^ip (p2 + m2)2 - ^^ ip (p2 + ^2)2 d^^^^^^P) + ^^2n8.-2 • (81) 



C. Gap Equation 

Upon inserting our solutions to the variational equation into (|57| ) and expanding it in 
powers of m/T, we will obtain a thermodynamic potential ^o{fn) that depends on a single 
variational parameter m. This function will contain severe ultraviolet divergences. Our 
goal will be to construct a thermodynamic potential that is a finite function of m but has 
the same variational minimum. To accomplish this, we will add to f2o(^^) a quantity that 
includes a factor of {m'^ — (7)2, where G is the right side of the gap equation (0). We will 
therefore need this gap equation in the form of an explicit equation for m. 

To make the gap equation (|62D explicit, we must expand the sum-integrals on the right 
side in powers of m/T. We will keep all finite terms through order g^ and all divergent terms 
through order g"^. Our first step will be to expand them in powers of the self-energy functions 
^n,k{P) and an,k{p)- Each of the sum-integrals in (^) can be expressed as a sum of the 
corresponding massive sum-integral, sum-integrals over P involving the functions n„^fc(P), 
and integrals over p involving the functions (Tn,k{p)- The terms that contribute to the gap 
equation through order g"^ are 

(r4,-2{q) + (T4,o{q) , ^ n4,o(Q) + n4,i(Q)' 



tp P2 + U(P) -^'~ ^^'^~^^' ^L ""'^'{ah rJ)2 ^ + 
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^screen -^screen '-'\90h^ ) 



^V7^^^W|k + q|) 



q (g2 + m2)2 



^^ , (X4,„2(g) 



k=im 

4 



+ 0{g') , (83) 



/q(g2+m2)2JpQ(g2)2 

The resulting form of the gap equation, including all term through order g'^T'^, is 
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(84) 



^ 2 yqV(g2 + m2)2 (g2+^2)3^^ 

The m/T expansions of the massive sum-integrals Xi and Xscreen are given in ( p.6|) and 
( p.9| ). The only momentum integrals involving self-energy functions that contribute to the 
finite terms through order g^T"^ or to divergent terms through order g'^T'^ are 
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The sum-integrals involving self-energy functions are 
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Keeping terms of order e through order g^, finite terms through order g^, and divergent 
terms through order g'^, the gap equation ( pi] ) becomes 
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D. Thermodynamic Potential 

The 3-loop ^-derivable approximation to the thermodynamic potential fio[-D] is given in 
(^) as a functional of the self-energy n(P). The sum-integrals can be expanded in powers of 
771 /T using the methods of Appendix 0. The function 11 (P) can then be eliminated in favor of 
the mass variable 7n by using the solution to the variational equation in subsection 0.B. This 
reduces the thermodynamic potential to a function f2o(m) of a single variational parameter 
m. This function contains many ultraviolet divergences. Some of them are eliminated 
when f2o(m) is evaluated at the solution of the gap equation. Some can be eliminated 
by renormalization of the coupling constant. There are other ultraviolet divergences that 
cannot be canceled and represent unavoidable ambiguities in the ^-derivable approximation. 
We wish to construct a thermodynamic potential that is finite to as high an order in g as 
possible. 

Our first step in constructing a finite thermodynamic potential fl{m) will be to cancel 
the most severe divergences that are eliminated by evaluating Qo{m) at the solution to 
the gap equation. These are the divergences proportional to m^/e and ^fgrri^/e^ that were 
already encountered in the 2-loop ^-derivable approximation. They can be canceled without 
changing the free energy by adding the term 
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Since this is proportional to the square of the gap equation evaluated at P = 0, it vanishes 
at the solution to the gap equation. Using (|6^) to eliminate n(0) in favor of m, this term 
can be written 
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The thermodynamic potential fli{m) defined by adding (|91|) to (|57| ) is 
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We want to expand the thermodynamic potential (192|) in powers of m/T, keeping all 
finite terms through order g^ and all divergent terms through order g'^. Our first step 
will be to expand in powers of the self-energy functions n„jc(P) and (Jn,k{p)- Each of the 
sum-integrals in (|92D can be expressed as a sum of the corresponding massive sum-integral, 
sum-integrals over P involving the functions n„^fc(P), and integrals over p involving the 
functions <Jn,k{p)- For the tadpole sum-integral and for Onscreen, these expressions are given 
in (|S^ ) and (|53|). For the other two sum-integrals in (^), the terms that contribute to the 
thermodynamic potential through order g"^ are 
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The resulting expression for the thermodynamic potential, including all terms through order 



^^ is 
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The next step is to express the coefficient of each power of go as an expansion in powers of 
m/T. The expansions of the sum-integrals Xq, Xi, Xgcreen, and Xbaii are given in Appendix 0. 
The integrals involving the self-energy function cr4 _2(g) are given explicitly in (|78D-(pTD. The 
only integrals that contribute to the finite terms in ( ]95|) through order g^ or to divergent 
terms through order (^^ are ( ]85| ) and 
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The complete m/T expansion of the thermodynamic potential, including all finite terms 
through order g^ and all divergent terms through order g''^, is 
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This expression has ultraviolet divergences of order g"^ and higher. We have two ways to 
eliminate these divergences. One is to use (^ to eliminate the bare coupling constant ao in 
favor of the renormalized coupling constant a{fi). The other is to add to the thermodynamic 
potential i7(m) something proportional to the square of the gap equation. In particular, we 
can add a term such as (0), with a coefficient proportional to any power of {gQfi~'^Y. We 
can define a thermodynamic potential 1^2 (m) in which the divergences are postponed to 
order g^ as follows: 
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where G is the right side of the gap equation (|89|). The factor of /i^*^ in the term that 
is added guarantees that 1^2 (^^i) will be independent of fi. The price that must be paid 
for the independence of fi is the introduction of another arbitrary scale A. After making 
the substitution (^ and truncating after terms of 5*^ order in rh or g, we obtain a finite 
thermodynamic potential 
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where i = log(A/47rT). This is our final result for the 3- loop ^-derivable thermodynamic 
potential. The divergent terms of order g^ that were eliminated by the truncation are 
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FIG. 4. <l>-derivable approximation to the screening mass divided by the leading order solution 
TO-Lo — 5'^(27rT)T^/24 as a function of g{2iTT) for (a) A = 2itT and (b) K = m. Solid lines are 
solutions to the gap equation truncated at orders (/", n = 2, 3, 4, 5 with ^ = 2itT. Dashed lines are 
solutions to the gap equation truncated at order g^ for fi = ttT and fi = AttT. 
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Note that the double poles in e do not cancel if we use the leading order gap equation 
m^ = a/6. Thus they cannot be eliminated by adding another term to the thermodynamic 
potential proportional to (m^ — G)"^. These ultraviolet divergences arise from having trun- 
cated the interaction functional $ at 3 loops. To eliminate them, one would have to include 
4- loop terms in $. 

In the case of the 2-loop <l>-derivable approximation, we were able to eliminate the 
ultraviolet divergences to all orders by replacing the true coupling constant a by a new 
coupling constant a whose expansion in a had divergent coefficients but which still had a 
well-defined beta function. If we demand that the resulting beta function has the correct 
one- loop coefficient, the new coupling constant must have the form 
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where A is a constant. By an appropriate choice of this constant, one might be able to 
ehminate one of the terms of order g^ with a single pole in e. But since we cannot eliminate 
the double poles, there is nothing to be gained by introducing such a coupling constant. 

E. Screening Mass 



The gap equation obtained by varying our finite thermodynamic potential (|99D with 
respect to m^ while holding /i and A fixed is 
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The solution to this equation is the 3-loop ^-derivable approximation to the screening mass 
rUs. If we solve the equation iteratively in powers of (a/6)^/^, it agrees with the weak-coupling 



expansion (M) for the screening mass. The screening mass depends on the arbitrary scale 
A through the logarithms i. It also depends on the renormalization scale /i through a{fi) 
and L. However since it agrees with the weak-coupling expansion through order g^, all 
dependence on these arbitrary scales must be of order g^ or higher. 

We can obtain a series of successive approximations to the screening mass by truncating 
the gap equation ( p.02| ) after terms of n*^ order in g and rh. In Fig. 4a, we show the solutions 
for m'^ divided by the leading order result ttlIq = (7^(27rT)T^/24 for yU = A = 2'kT as a 
function of g{2iTT). The solid lines correspond to truncations after orders n = 2, 3, 4, 5. For 
n = 4 and 5, there is no solution beyond a critical value of the coupling constant: g = 2.591 
for n = A and g = 2.768 for n = 5. For larger values of g, the thermodynamic potential 
has a runaway minimum m ^ that comes from the a^logm term in (p9|). A similar 



phenomenum occurs in screened perturbation theory at three loops |2T|: if the screening 
mass is used as the variational parameter, the solution to the gap equation terminates at 
nearly the same critical value of g. For g below the critical value, the screening mass seems 
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FIG. 5. <I>-derivable approximation to the pressure divided by that of the ideal gas as a function 
of g{2TTT) for (a) A = 2ttT and (b) A = m. Solid lines correspond to the truncation of the 
thermodynamic potential at order g^, n = 2,3,4,5 for fj, = 2'kT. Dashed lines correspond to the 
truncation at order g^ for /i = vrT and ^ = 47rT. 

to be converging nicely. It is also fairly insensitive to the renormalization scale, as shown in 
Fig. 4a. The solution to the gap equation truncated at 5*"^ order is shown as dashed lines 
for the two choices /x = vrT and ^ = AttT. 

The simplest way to get a solution for the screening mass that extends beyond the 
critical value of g is to choose the arbitrary scale A to be proportional to m. This changes 
the gap equation to 
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When truncated at 2"^^ and 3^^^^ order in g and rh, the gap equations ( [I02| ) and ( |103| ) are 
identical. The solution to the gap equation ( |103D truncated at n*^ order and with /x = 2nT 
and A = m are shown as solid lines in Fig. 4b. The solution for n = 4 terminates at 
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FIG. 6. <I>-derivable approximation to the entropy divided by that of the ideal gas as a function 
of g{2TTT) for (a) A = 2'kT and (b) A = m. Sohd hnes correspond to the truncation of the entropy 
at order g"^, n = 2,3,4,5 for /i = 2'kT. Dashed hnes correspond to the truncation at order g^ for 
/i = vrT and /i = 4ttT. 

g = 3.551 for fi = 2'kT . The solution for n = b continues well beyond g{2'KT) = 4. The 
dashed lines are the 5*^^ order truncation with A = m for two choices of the renormalization 
scale: /i = vrT and /i = AttT. 

F. Thermodynamic Functions 

The pressure, given by (||), is obtained by inserting the solution to the gap equation ( |102| ) 
into the thermodynamic potential (pQ]). The pressure truncated at n^^ order in g and rh is 
obtained by using the n^^ order truncations of ( |102| ) and (|99|). In Fig. 5a, we show the n^^ 
order truncations of the pressure divided by that of the ideal gas for A = 2ttT and /i = 27rT. 
We show also the change in the 5**^ order truncation from varying the renormalization scale 
fi. Since the solutions to the 4**^ order and 5*^^ order gap equations cannot be continued 
beyond the critical values g{2'KT) = 2.591 and 2.768 respectively, the corresponding lines for 
the pressure cannot be extended past these values. The variations of the pressure are much 
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smaller than one might expect from the variations of the screening mass in Fig. 4a. We see 
an obvious improvement in the convergence compared to the weak-coupling expansion. At 
gth Qj-fjgj.^ gi^^d up to the critical values of g{2TTT) starting around 2.5, the variations from the 
renormalization scale are much smaller than those of the weak-coupling expansion. Fig. 5b 
shows the pressure truncated at the same orders but with the scale A = m. The pressure 
is obtained by inserting the solution to the gap equation (|103|) into the thermodynamic 
potential (|99|) . The line for n = 5 continues well beyond g{2iTT) = 4, but for values of 
g{2TTT) around 2.5 the convergence between n = 3, 4 and 5 is not as good as in the previous 
case of A = 2'kT. Up to g{27fT) = 4, the changes in the 5*^^ order truncation from varying the 
renormalization scale /z appears to be much smaller than in the weak-coupling expansion. 

The diagrammatic entropy is given by Eq. (P), with /i and a(yu) held fixed. We have to 
consider two cases giving different expressions of the entropy, the case with A = c(27cT) and 
the case with A = c m. In the first case, upon differentiating the thermodynamic potential 
(p9|) with respect to T, the 5*^^ order expression for the entropy is 
1/5 \ 1 r 
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where iSidcai = (27r^/45)T^, and m is the solution of the gap equation ( |102| ). Fig. 6a displays 

the good convergence of the ra*^ order truncation of the entropy for fi = 2ttT and A = 27rT, 

and the small variations from the renormalization scale /x. The 5**^ order expression for the 

entropy for the choice A = c m is given by 
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with rh the solution to the gap equation ( |103| ) and i = log(m/27rT). Fig. 6b shows the 
^th Qj^jgj^ truncations of the entropy for A = m. We conclude again, as in the case of the 
pressure, that around the critical values of g{27fT) the convergence is not as good as in 
the case of A = 2nT. However the changes of the 5*^ order truncation with respect to the 
variations of the renormalization scale fi are still relatively small. 

VI. CONCLUSIONS 

We have developed a systematic method for solving ^-derivable approximations for the 
massless 0^ field theory. The method involves expanding sum-integrals systematically in 
powers of m/T, where m is the screening mass. When expanded in powers of m/T, the 
variational equations for the self-energy n(P) have a recursive structure that allows then 
to be solved in terms of the variable m. This reduces the problem to one that involves 
the single variational parameter m. Upon truncating the expansion in m/T, we obtain a 
gap equation for m that could easily be solved numerically were it not for the presence of 
ultraviolet divergences. Some of the ultraviolet divergences vanish at the variational point, 
and can be eliminated by a redefinition of the thermodynamic potential. However there are 
other ultraviolet divergences that do not vanish at the variational point. They are generated 
by the truncation error in defining the interaction functional $[H]. They limit the accuracy 
to which the ^-derivable approximation can be solved. 

In the case of the 2-loop ^-derivable approximation, the unavoidable ultraviolet diver- 
gences appear at order g'^. They can be eliminated to all orders by a trick that involves 
introducing a new coupling constant g that runs with the wrong beta function. However 
this trick does not generalize to higher orders in the loop expansion. If one expresses the 
thermodynamic potential in terms of the true coupling constant g, the best one can do in 
constructing a finite thermodynamic potential is to truncate after terms of third order in 
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m/T. 

In the 3-loop ^-derivable approximation, the unavoidable uhraviolet divergences appear 
at order g^. We constructed a finite thermodynamic potential by adding terms proportional 
to the square of the gap equation to eliminate the ultraviolet divergences that vanish at the 
variational point and then truncating after terms of S*'^ order in g and m/T. Truncations 
of the finite thermodynamic potential at n^^ order in g and m/T ., where n = 2, 3, 4, 5, 
define a sequence of successive variational approximations to the thermodynamics. The 
truncation after terms of 3'''^ order in g and m/T reproduces the finite thermodynamic 
potential defined by the 2-loop ^-derivable approximation. Our solution of the 3-loop <l>- 
derivable approximation allows the 4**^ order and 5^^ order approximations to be constructed. 
The 6^^ order and 7**^ order approximations could be constructed by applying our method 
to the 4-loop ^-derivable approximation. The only obstacle is evaluating the scalar sum- 
integrals that would arise. The most difficult of the sum-integrals are the same ones that 
would need to be evaluated to extend the weak-coupling expansion to order g^. 

We calculated numerically the thermodynamic functions defined by the finite thermo- 
dynamic potential truncated at n^^ order in g and m/T for n = 2,3, 4, and 5. One problem 
is that in the cases n = A and 5, the solution to the gap equation for m cannot be extended 
beyond a critical value of g that is near g = 2.5 if we choose the renormalization scale 
fi = 2ttT . For values of g below the critical value, the successive approximations seem to be 
converging and are rather insensitive to variations in the renormalization scale. 

Our construction of a finite thermodynamic potential involved a prescription for remov- 
ing from the bare thermodynamic potential VLq ultraviolet divergences that vanish at the 
variational point. Our prescription was to add the two terms in ([91| ) and (p8|) that are both 



proportional to the square of the gap equation. We had no compelling motivation for this 
prescription beyond its success in the 2-loop ^-derivable approximation. This prescription 
introduces an ambiguity into the finite thermodynamic potential (pQ]). We can feel free to 
add additional terms proportional to (m^ — G)^, where G is the right side of the finite gap 
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equation given in ( [L02| ) or (|103|) . If we avoid large coefficients for these additional terms, 
they do not have a large effect numerically for g less than about 2.5. However the ambiguity 
emphasizes that our prescription was rather ad hoc. It would be desirable to have a much 
more compelling prescription for eliminating the ultraviolet divergences that vanish at the 
variational point. 

Our strategy for solving ^-derivable approximations should also be applicable to more 
complicated field theories, such as gauge theories. In the case of QCD, Blaizot, lancu 
and Rebhan |jT^ have solved the 2-loop ^-derivable equations for the self-energy tensor 



n^'^ to leading order in Qs in terms of a Debye screening mass parameter m^). They did 
not determine m£, variationally, but simply used the asymptotic weak-coupling expression 
mu = QgT. The resulting approximation reproduces the thermodynamic functions through 
order g"^. They called it a "leading order" (LO) HTL approximation. The application of our 
method to gauge theories would involve expanding sum-integrals systematically in powers 
of Qs and rriij/T. If applied to the 2-loop <l>-derivable approximation for QCD, it should 



reproduce the leading order solution for I[^'^{P) in Ref. |]n|, but it should give a gap equation 
for niD that allows it to be determined variationally. By truncating after terms of 3"^*^ order 
in Qs and mo/T, one should be able to improve the accuracy of the thermodynamic functions 
to order g^^. At higher orders in Qs and mo/T, the 2-loop ^-derivable approximation will 
be plagued by ultraviolet divergences associated with the truncation of $ as well as by 
gauge dependence associated with the use of the gauge-dependent gluon propagator as a 
variational function. To improve the accuracy while avoiding these problems, one would 
have to apply our method to the 3-loop ^-derivable approximation to QCD. This would 
allow the accuracy to be improved to orders g^ and g^. Any further improvement in the 
accuracy would be very difficult. One problem is that there are new sum-integrals that enter 
at order g^ that have never been evaluated. A more serious problem is that the 3-loop free 
energy is sensitive to the nonperturbative momentum scale g'^T associated with magnetic 
screening. This will probably cause our method for reducing the ^-derivable approximation 
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to a single-parameter variational problem to breakdown. 



Blaizot, lancu, and Rebhan [T^] have introduced an improvement on their LO HTL 
approximation that reproduces the thermodynamic functions to order g^. They call it a 
"next-to-leading order" (NLO) HTL approximation. It reduces essentially to the substitu- 
tion 

^l — ,!_\ , (106) 

in their LO HTL approximation. The simpler prescription mfj -^ glT'^{l — Sqs/tt) gives the 
same parametric improvement in the accuracy, but is a disaster phenomenologically, giving a 
negative screening mass unless a^ < 0.087. The ad hoc prescription ( |106| ) gives phenomeno- 
logically acceptable results that are not very different from the LO HTL approximation of 
Ref . Jl^ . We argue that such a prescription is unlikely to arise in any systematically improv- 



able approximation to the thermodynamics. Our solution to the ^-derivable approximation 

to the massless 0^ field theory suggests an alternative prescription. The truncation of a gap 

equation after terms that are S''*^ order in Qs and m£,/T as in ( p6D would give 

9 .\i/2 3 
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In summary, we have developed a new method for solving the 3-loop ^-derivable approx- 
imation for the massless 0^ field theory. This method can be used to construct a sequence of 
systematically improvable approximations to the thermodynamic functions. This sequence 
of approximations seems to be stable and converging for g{2'n'T) less than 2.5, and also 
insensitive to the choice of renormalization scale. If this method can be adapted to gauge 
theories, it may provide a practical solution to the longstanding problem of large radiative 
corrections to the thermodynamic functions. 
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APPENDIX A: EXPANSIONS IN M/T 

In this Appendix, we carry out the expansions of sum-integrals in m/T. We first illus- 
trate the method by applying it to one-loop sum-integrals for a particle of mass m. We then 
apply it to sum-integrals with a momentum dependent self-energy n(-P). 

a. Simple examples 

To illustrate the m/T expansion in the simplest context, we carry it out for the one-loop 
sum-integrals that appear in the 2-loop <l>-derivable approximation. The sum-integrals over 
P = {27mT, p) in ( pi]) and (^) include sums over all Matsubara frequencies n and integrals 
over all momenta p. There are two important momentum scales: the hard scale 27rT and the 
soft scale m. The hard region for the momentum P = (27rnT, p) includes n 7^ for all p and 
also n = with p of order T. The soft region is n = and p of order m. We will evaluate 
the contributions to the sum-integrals from these two momentum regions separately. 

The tadpole sum-integral in p^ ) can be separated into contributions from hard and soft 
regions, which we denote by the superscript (h) and (s): 



1 rW 1 ^(s) 1 



fp P'^ + m? jp P'^ + im? ' Jp P2 + m^ 
In the hard region, P^ is of order T^, so we can expand in powers of m? / P'^: 



(A.l) 



Tp P^ + m^ Tp\P^ (P2)2^£-2 {P^Y+^I' ^ ' ' 

The contribution from the soft region can be written 

(^) 1 /■ 1 

T -r—3- (A-3) 



Tp P"^ + m? Jp p^ + rri 

One could accomplish the separation of the hard and soft momentum regions by introducing 
a separation scale A satisfying m ^ A ^ 2ttT. The dependence on A would then cancel 
between ( |A.2[ ) and ( |A.3D . However, if we use dimensional regularization, we can get the 
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correct answer more easily by taking the limit A — i> in ( |A.2| ) and A — i> oo in ( [A.3|) . The 



complete expansion of the sum-integral in powers of m/T is then just the sum of (|A.2| ) and 



1 



1 /• 1 °° /■ 1 



fpP^ + m? fpP^ ' % (P2)2 ' ^2 ^P (P^ 

where /i is the momentum integral on the right side of ( [A .31 ). 



The logarithmic sum-integral in (^) can also be separated into contributions from hard 
and soft regions, as in ( |A.1|) . In the hard contribution, we can expand in powers of rn? /P"^: 

f^ log(P^ + ™^) = f^ ^.og P= + ^ - -^^ + g '-^j^) . (A.5) 

The soft contribution is 

i ' log(p2 + m^) =T f log(p2 + m^) . (A.6) 

The complete expansion of the sum-integral in powers of m/T is the sum of ([A .51) and ( |A.6| ): 



^ l0g(p2 + ^2^ = ^ log P2 + ^2^ i_ ^ y (_ J. ) 



where — /g is the momentum integral on the right side of ( |A.6|) . 

b. Tadpole sum-integral 

The simplest sum-integral with a momentum-dependent self-energy is the tadpole sum- 
integral, which appears in the expression ( |62D for the screening mass. The tadpole sum- 
integral can be separated into contributions from hard and soft regions, which we denote by 
a superscript (h) or (s): 

r 1 _rW 1 .(s) 1 

Tq q2 + n(Q) ~ ¥q q2 + n(Q) + jfe Q2 + n(Q) ' ^ ' ^ 

In the hard region, Q^ is of order T^ while 11 is of order m^, so we can expand in powers of 
n((5^). Inserting the expansion (^), we obtain 
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{goi^-rt "^•°^^,^+,^^^^^ + 0(/) . (A.9) 



Q g2 + n(g) ^q VQ' W')' W')' 

O + n 

The first term on the right side is proportional to T^"^*^. Taking m to be of order g, the 
others are suppressed by g"^, g^, g^, and g^, respectively, in the limit e ^ 0. The contribution 
from the soft region can be written 



'Iq Q2 + n(g) ^ J^q^ + m' + a{q)- ^^'^^^ 

In this region, q^ is of order m?, which is the first term in the expansion (|66|) of 11(0, p). 
Expanding in powers of the higher terms, we have 



■J<\\{q^ + rn?Y {q^ + m?Y ) 



The first term on the right side is proportional to in^~'^'^T. In the limit e — > 0, it is of order 
g when m is of order gT . In this limit the other terms are of order g^ , g^ , g^ , and g^, 
respectively. The error estimate 0{g'^) on the sum-integrals ( |A.9D and ( |A.11D refers to their 
orders in g in the limit e — > 0. 

The complete expression for the tadpole sum-integral is the sum of the hard and soft 
contributions in ( [A.9|) and ( |A.11|) . The integrals and sum-integrals that do not involve the 



self energy functions Yin,m{Q) and <Jn,m{q) are given in Appendices ^ and 0. 



c. Sunset sum-integral for hard P 



The sunset sum-integral ( |60D depends on the external momentum P, and the expansion 
in powers of m/T is different for hard P and soft P. We first consider the case of hard 
P. There are 3 different sum-integration regions that must be considered. The momenta Q 
and R can both be hard, one can be hard and the other soft, or they both can be soft. We 
denote these regions by superscripts (hh), (hs), and (ss), respectively. 
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In the (hh) region, we can expand in powers of the self-energy. Only the m'^ term in the 
self-energy (^) contributes up to order g^T^: 

We next consider the (hs) region for hard P. Taking into account that any of the three 
momenta Q, R, and S = —{P + Q + R) can be soft, the sum-integral can be written 

^sun (P) = 32^ / J g2 + n(Q) r^ + m^ + a(r) (P + Q)2 + 2(p + q) -r + r^ + U{S) ' ^^'^^^ 

where S = —{P + Q) — (0,r). The m/T expansion is obtained by expanding the sum- 
integrand around l/[Q'^{r'^ + m?){P + QY]. The first few terms are 

With dimensional regularization, the integral Jq vanishes. The factor of 1/d multiplying one 

of the sum-integrals comes from averaging over the angles of r in (i dimensions: [(p+q) -r]^ — >■ 

(p + q)^r^/(i. The change of variables Q -^ —P — Q then puts the sum-integral over Q into 

the form shown. 

Finally, we consider the (ss) region for hard P. Taking into account that the hard 

momentum can flow through any of the three propagators, the sum-integral can be written 

c-(ss)/p\ ^ 07.2 /■ \ \ \ 

^""^ ^ yqrg2 + m2 + (T(g)r2 + m2 + a(r)P2 + 2p-(q + r) + (q + r)2 + n(5)' 

(A.15) 

where S = —P — (0, q + r). The m/T expansion is obtained by expanding the sum-integral 
around l/[(q'^ + m'^){r'^ + m^)P^]. The leading term is of order m^ and all higher terms are 
of order g^ or higher. Thus the {ss) term in the sum-integral can be written 

'^iuXP) = 3T'lf^^+0{g% (A.16) 
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The complete expression for the sunset sum-integral for hard momentum P, including all 
terms up to errors of order g^, is the sum of ( |A.12| ), ( |A.14| ), and (|A.16| ). 



d. Sunset sum-integral for soft P 

We next carry out the m/T expansion for the sunset sum-integral with soft external 
momentum P = (0, p), where p is of order m. Again, there are three momentum regions to 
consider: (hh), (hs), and (ss). 

We first consider the (hh) region. The sum-integral can be written 
r I I 1 

^sun (0, p) = J^^ Q2 + n(Q) i?2 + u{R) (Q + ny + 2(q + r) ■ p + p^ + u{s) ' ^^-^^^ 

where S = ~{Q + R) — (0,p). The m/T expansion is obtained by expanding the sum- 
integrand around 1/[Q'^R'^{Q + i?)^]. The first few terms are 

^iunHO, P) = f^^ g2^2(g + ^)2 - iP' + 3"^')fQ^ (Q2)2^2(g + ^)2 



r TQRiQ^yR^iQ + R)^ 
The first sum-integral on the right side vanishes. 

We next consider the {hs) region. There must be one propagator with soft momentum. 
Taking into account that it can be any one of the three propagators, the sum-integral can 
be written 

^"'^^ '^^ yr^g2 + n(g)r2 + m2 + a(r)g2 + 2q.(p + r) + (p + r)2 + n(5) ' 

(A.19) 

where S = — Q — (0, p + r). The m/T expansion is obtained by expanding the sum-integrand 
around 1 / [{Q'^)'^ {r'^ + m"^)]. The first few terms are 

3S:?(0, P) = STI^^ j^_ - ST [/„ + (p= + ™=)A] i (^ 






r^ + m^yTQ (Q2)2 



+ T^ 



Jo + ip' - m')h\ f^ ^ + Oig') . (A.20) 



45 



The integrals Iq vanish in dimensional regularization. 

Finally, we consider the (ss) region. The sum-integral in this region can be written 

1 1 1 






(A.21) 



/qr q^ + m? + (7{q) r^ + m? + a{r) s^ + m? + a{s) 

where s = — (p + q + r). The m/T expansion is obtained by expanding in powers of cr(g), 
cr(r), and cr(s). The first few terms in the expansion are 



S(-)(0, p) = T'hunip) - 3igof^-rT' f ^/2'^t2 -^bub(|p + q|) + 0{g') , (A.22) 



where Isunip) and /bub(p) are the sunset and bubble integrals defined in ( |(J.1(J| ) and ( |C9D . 

The complete expression for the sunset sum-integral at soft momentum P = (0,p), 
including all terms up to errors of order g"^, is the sum of (|A.18| ), ( |A.20| ), and ( |A.22D . 
Evaluated at p = im, it reduces to 



^sun(0,p) 



T=/.„.,(™) + 3TI,f^ ^ - f /,™^T^^ j|jj 



-2m 



g2 + (2/(i)q2 



QR {Q^)m^{Q + Ry 

Hgofi-r \t' I T^^f^/bubdp + q|) 
Jq [q'^ + ni'^Y 



+T 



p=im 



C^4,-2('") 



{r^ + m?fTQ (Q2) 



+ 0{g^) 



(A.23) 



e. Sum-integrals in the thermodynamic potential 



We now apply the m/T expansions to the sum-integrals that appear in the thermody- 
namic potential (|57D . The method should be clear from the example of the sunset diagram, 
so we only give the results for each momentum region. 



For the logarithmic sum- integral in (57), the contributions from the hard and soft regions 



are 



f^\og{P^ + Ii{P)) = f^ 



TTl 



logP^ + -;r 



m 



+ 



m 



p2 2(P2)2 3(p2 
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+ 



(gof^'Tf^ 



n4,o(P) + n4,i(P) + n4,2(P) + n4,3(P) 



-m 



p2 

2n4,o(P) + n4,i(P) 



+ i9ofi-'ff^ 



(P2)2 



p2 



+ 0{g' 



r-(s) 



(^4-2{p) +Cr4,o(p) 



^^' iog(P' + n(P)) = T{-Q + {gof^-rTJ^ 

' 0-8,-4 (P) ^1,-2 (P) 



p2 _|_ j^2 



+ ((7o/i-^)'T / 



+ (^(/ 



(A.24) 



(A.25) 



'P \P^ + "^^ 2(p2 -|- 77^2)2^ 

For the other one-loop sum-integral in the thermodynamic potential, the hard and soft 
contributions are 



w n(P) 



p p2 + n(P) 



m 



m 



P2 (P 



2^2 



+ 



m 



2^3 



■(fl'o/i 



-e\4 



(P 

n4,o(P) + n4,i(P) + n4,2(P) + n4,3(P) 



p2 

2n4,o(P) + n4,i(P) 



+ i9o^i'l'f^ 



2"^ (•p2p 

n8,-2(P) + n8,-i(P) 



p2 

2(^4,~2{p) + (TiM 



+ 0{g') , 



W n(P) r 

P p^ + n(P) = ^ ^^^ + (^0/^ )Ty^^ 



(p2 -|- m'^y 

p' 



2 a8,-4(p) „2 ^1,-2 (P) 



c?(^' 



(A.26) 



(A.27) 



The two-loop sum-integral in (0) is the square of the tadpole sum-integral. It is mul- 
tiplied by Qq, so to get the free energy through order (7^, we need the square of the tadpole 
integral to order g^. This is obtained by squaring the sum of ( |A.9|) and (|A.ll ), and keeping 
all terms to the desired order. 

The basketball sum-integral in (|57|) involves a triple sum-integral, so there are 4 momen- 
tum regions: (sss), (hss), (hhs), and (hhh). The contribution from each of these regions 



IS 



^ 



{hhh) 

ball 



Aw? 



PQR \P^Q^R\P + Q + RY {p^^qm^iP + Q + Ry 



+ 0{g' 



(A.28) 



<y 



(hhs) 
ball 



4TJii 



QRQm\Q + Ry 



AT (jo + 2m^h) i 



QRiQ-^YR^Q + Rf 
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16 / 2t\J< q^ 



C^{hss) 
^ball 



^ball 



rf V " '^ TqR (Q2)3i?2(g + ^)2 

QT^llf^ (^ + ^^3') , (A.30) 

T^/baii - K9.^^-rT' I ^t'^^Ti^^ + ^(^') ' (A-31) 



where Jbaii is the basketball integral given in ( |C.13| ). In ( [A.29| ), the first sum-integral and 
the integral Iq vanish with dimensional regularization. 

APPENDIX B: MASSLESS SUM-INTEGRALS 

In this Appendix, we list the massless sum-integrals that appear in the calculation of 
the 3-loop ^-derivable free energy through 7**^ order in g and m/T. Analytic expressions 
are given for all the sum-integrals with the exception of one 3-loop sum-integral, for which 
the constant term has not been evaluated. 

In the imaginary-time formalism for thermal field theory, a boson has Euclidean 4- 
momentum P = {po,p), with P^ = pp + P^- The Euclidean energy po has discrete values: 
Po = 27mT, where n is an integer. Loop diagrams involve sums over po and integrals over p. 
With dimensional regularization, the integral is generalized to (i = 3 — 2e spacial dimensions. 
We define the dimensionally regularized sum-integral by 

where c? = 3 — 2e is the dimension of space and fi is an arbitrary momentum scale. The factor 
{e"^ /Airy is introduced so that, after minimal subtraction of the poles in e due to ultraviolet 
divergences, /i coincides with the renormalization scale of the MS renormalization scheme. 
The basic one-loop sum-integral is 

1 C(2ri - 3 + 2e) r(r2 - I + e) 



P (P2)n 87r2 r(i)r(n) 
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(eV)'(27rT)^-2"-2e^ (32) 



where 7 is Euler's constant and ({z) is Riemann's zeta function. Most of the exphcit one- loop 
sum-integrals required in our calculations can be derived from ( |B.2|) : 



logP' 

1 
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12 V2 



[l + 0(e)] 



1+ ( 2 + 2^^ 1 e + 0(6^) 



C(-i^ 



P(P2)2 (47r)2 1^2^ 

1 1 



2e 



■K 



- + 27+ — -471 e + 0(e^ 



[2C(3)+0(6)] 



(B.3) 
(B.4) 

(B.5) 
(B.6) 



fp (P2)3 (4vr)4r2 

where fi = /i/(27rT). The number 71 is the first Stieltjes gamma constant defined by the 
equation 



C(l + z) = - + 7 - 7i^ + 0{z^). 



We also need one other one-loop sum-integral: 

1 



Pi 



P (P2)4 (47r)4T2 



[C(3) + 0(e)] 



The two-loop sum-integrals that are needed are 
1 



PQP^Q^iP + qy 
1 







fX 



4e 



PQ(P2)2g2(p + g)2 2(47r)4 V2, 



■^ + (1+47)- 



IT 



+ ( 3 + y + 47 + 472 - 871 J + 0(e) 



pI 



1 I ji 



4e 



PQ(P2)3Q2(p + g)2 8(47r)4 ^2 



1 /9 \ 1 



/35 7r2 \ 

+ y + Y + I87 + 47' - 871 + 0(e) 



The sum- integrals (p.lOp and ( p.ll ) can be evaluated using the methods of Ref. 
will need two specific hnear combinations of these sum-integrals: 



P2 + (2/rf)p2 ^ 3 f^' 
PQ{P^fQ\P + QY ~4(47r)4 I 2^ 



4e 



1 [^ ^ \^ 



(B.7) 



(B.8) 



(B.9) 



(B.IO) 



(B.ll) 



1. We 



49 



/^89 7r2 10 o \ ,^, , 



P2 - (4/t/)p2 



/i 



4e 



PQ {P^fQ^P + Q)2 4(47r) 



^+(y + 47)+0(6) 



(B.12) 
(B.13) 



The three-loop sum-integrals that are required are 

6e 



1 



rpi 



PQRP^qm^S^ 24(47r)2 1^2^ 

1 ^ _T^ /M 

pqr{PYQ^R^S^ 8(47r)4 1^2/ 



91 C'(-l 



1 /89 , 
^+ y + 4, + 2 



C(-3), 
C(-1)V 

c(-i);e. 



0(e) 



(B.14) 
(B.15) 



where S" = P + Q + -R- The niassless basketball sum-integral (P.14| ) was first evaluated 
analytically in Ref. [^]. The sum-integral (p.l5| ) is more difficult. One way to evaluate it is 



to use the result for the massive basketball sum-integral in Ref. pO[- ^^ ^^at sum- integral is 
expanded in powers of m/T as in (p.ll|) , the sum- integral ( |B.15| ) appears at order {m/T)'^. 
One can deduce the pole terms from the analytic expressions for the pole terms in the 
massive basketball sum-integral in Ref. [^. We have not evaluated the constant term in 
this sum-integral. 

The transcendental constants that appear in these sum-integrals are 



7 =0.57722, 

7^ = -0.072816; 

C(3) = 1.20206, 

C'(-l)/C(-l) = 1.98505, 

C'(-3)/C(-3) = 0.64543. 



(B.16) 
(B.17) 
(B.18) 
(B.19) 
(B.20) 



APPENDIX C: INTEGRALS 

In this Appendix, we list the integrals that appear in the calculation of the 3-loop $- 
derivable free energy through 7**^ order in g and m/T. The integrals that contribute through 
gth qy(Iqy g^]^^ those that contribute to the pole terms at 6**^ and 7*^^ order are evaluated 
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analytically. Some of the integrals that contribute to finite terms at 6^^ and 7**^ order are 
evaluated numerically. 

In a Euclidean field theory in 3 space dimensions, loop diagrams involve integrals over 
3-momenta. With dimensional regularization, the integral is generalized to d = 3 — 2e spacial 
dimensions. We define the measure for the dimensionally regularized integral by 



Air 



(27r) 



3-2e 



(C.l) 



If renormalization is accomplished by the minimal subtraction of poles in e, then /x is the 
renormalization scale in the MS scheme. 

The basic one-loop integral for a boson with mass m is 

^ (C.2) 



p (p2 + m^)"- 



Differentiating with respect to n and then setting n = 0, we obtain the integral 



I^= f log(p2 + m^) . 

JJ3 



(C.3) 



The analytic expression for the integral ( \C'2\ ) is 

1 r(n-^ + 



^vr r(i)rH 



(eV)' 



m 



3-2n-2e 



(C.4) 



The explicit one-loop integrals required in our calculations can be derived from (|C.4|) : 



/o = 0, 

/' = _!!^ fJL] 

" ° QiT \2mJ 

^ 47r \2m) 



2e 



2e 



l + ^e + 0(e^) 



TT 



1 



(-) 

8nm \2mJ 



2e 



l + 2e+ \A+ — \e^ + 0{e^) 



l + 0(e' 



(C.5) 
(C.6) 

(C.7) 
(C.8) 



The multi-loop integrals that are needed can be conveniently expressed in terms of the 
bubble and sunset integrals defined by 
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hun{p) 



q g2 + m^ (p + q)2 + m'^ 



qr g^ + m'^ r^ + m^ (p + q + r)^ + m^ 



(C.9) 
(CIO) 



We need the analytic continuation of the sunset integral to p = im. It is evaluated in Ref. 
through order e°, but we also need the order e term: 



Isnniifn) 



^ ('_^V'[l + (6-81og2) 



4(47r)2 \2mJ 
' 111'' 



+ 24 



6 



27 + 872 - 48 log 2 + 8 log^ 2 e + O (e^ 



(C.ll) 



The three-loop basketball integral can be expressed as 

/sun(p) 



'ball 



1^ _L m-^ 



p p^ + m 



(C.12) 



This integral is also evaluated in Ref. 0: 



m 



'ball 



(47r)3 V2m/ 



6e r 



(8-41og2) + 0(e) 



(C.13) 



The other three-loop integral is 

/sun(p) - Isun{ifn) 



[(1- log 2) + 0(e)]. 



(C.14) 



/p (p^ + TTi^)^ 2(47r)^?71 

This can be obtained by differentiating ( |C.13| ) with respect to m^ and also using ( |C.8|) and 
( p.ll| ). The integral (|C.14| ) is convergent for d = 3, so it can also be calculated directly 
without using dimensional regularization. In the limit e — > 0, the function in the numerator 



IS 



Isnnip) - Isnn{im) 



(477) 



3m p 1 p + 9m 

atan \ — log — 

p 3m 2 GAm.-^ 



(C.15) 



The 4-loop integral that is needed is 



q (g^ + m- 



2A2 



-/bub(|p + q|) 



p=im 



{A-nYm'' 



[-0.022048 + 0(e)] . (C.16) 



This integral can be evaluated directly in c? = 3 dimensions. The bubble integral is 
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4ub(|p + q|) — ^ -r\ — \ — ratan^^ . (C.17) 

47r|p + q| 2m 

After averaging over the angles of q, we can set p = im. Inserting the resulting function 
into ( p.l6| ) along with ( |C.15[ ), we obtain a one- dimensional integral that can be evaluated 



numerically. 

Finally there are several 5-loop integrals that are needed. Two of them are very simple: 

r [hunip]- Isuninn)]' ^ l^ 2519 + 0(e)] , (C.18) 

Jp {p^ + m?Y (47r)6m ^ ^ '^ ^ ' 

I [hnn{p) - hnJM)f ^ 1_^ [0.012004 + 0(e)] . (C.19) 

They can be evaluated directly in c? = 3 dimensions by inserting the expression (|C.15| ) for 



the function in the numerator. The other 5-loop integral is 



i ^=Pw^^^-<^' - Wf^ i«-^**"^ ^ '^w' ■ '^-^"^ 



In the limit e — ;> 0, the last factor in the integrand reduces to 



;hun{p) 



dm? (47r) 



3 p 

-atan- 



2mp 3m 



(C.21) 



The integral in ( P.20|) can be evaluated directly in (i = 3 dimensions by inserting the 
expressions ( |C.15|) and ( |C.21|) for the functions in the numerator. 



APPENDIX D: MASSIVE SUM-INTEGRALS 

To calculate the m/T expansion of the thermodynamic potential, we need to expand 
several massive sum-integrals in powers of m/T . The m/T expansions can be obtained 
from those in Appendix ^ simply by omitting terms containing the self-energy functions 
^n,k{P) and (Jn,k{p)- The remaining terms in the m/T expansion involve the sum- integrals 
in Appendix |B] and the integrals in Appendix 0. 

The two massive one-loop sum-integrals that we need are 
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x; = 2;iog(p2 + m2), 



Xi 



p2 _j_ ^ 



(D.l) 
(D.2) 



Their m/T expansions are given in ( [A .41) and (|A.7| ) 



J^ = T(-/^) + f^ 
Xi = T/i + 



m 



logP" + — 



+ 



p2 (p2)2 (p: 



p2 2(P2)2 3(p2)3 



+ C(m^ 



2 4 

+ 



2^3 



(D.3) 
(D.4) 



Inserting the sum-integrals from Appendix ^ and the integrals from Appendix |C|, we get 



° 45 I 2 



2e 



1 - ISm^ + 60m^ + — f'- + 27 ) m^ 



15 



C(3)m6 ^ Q^^^-^ 



(D.5) 



X, = 1t^ ( ^ 

12 \ 2 



2e 



1+ (2 + 2^^-#|e 



6fl + (-21og?fi + 2)elm 



C(-i 

3 ("i + 27") m^ + |c(3)m4 + 0{m 



(D.6) 



where m = m/{27rT) and /i = ^/{2ttT). The errors in the coefficients of m" are one order 
higher in e than the last term shown explicitly. 

The screening sum-integral for a particle of mass m is 

1 



T 



P=(0,p),p='im 



/qr (Q2 + m2)(i?2 + rn2)((p + Q + p)2 + ^2) 

The m/T expansion is obtained by omitting the terms in (|A.23| ) involving (T4 _2(p): 

1 



(D.7) 



T 



T'hUim) + 3r/i: 



2 



-2m¥ ^' + ^'/^^^^ 

^/?(Q2)3p2(g + p)2 



-—hm^T 



q- 



+ O(m^) . 



(D.8) 



The sum-integrals and integrals are given in Appendices ^ and |C[ Keeping finite terms 
through order m and divergent terms through order m^, we get 
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T 



1 rr.2 (A 

4(4vr)2 \2^ 



4e 



-6 

3 

~2 

The massive basketball sum-integral is 



+ (-41ogm + 6-81og2' 
- + (-21ogm + 2 + 27) 



m 



1 /5 \ 1 

+ ( - + 47 



TTl 



(D.9) 



J> 



ball 



1 



(D.IO) 



where S = P + Q + R. The ?7i/T expansion of this sum-integral is obtained by dropping the 
terms involving cr4„2(p) in the sum of ( |A.28| )-( jA.31| ): 

/ 1 Am"^ \ 



X, 



ball 



PQR yp^qm^s^ {p^ygm^s^ 

Q2 + (2/d)q2 



To. 



+ Oim") . 



:d.ii) 



The sum-integrals and integrals are given in Appendices ^ and |C[ Keeping finite terms 
through order m and divergent terms through order m^, we get 

6e 



X, 



T-14 



/i 



ball 



24(47r) 



1 , 91 , oC(-l) gC(-3) 
6 15^ C(-l) C(-3) 

6 log m + 8 — 4 log 2 



-12 
-3 

+1 



4+ ?+4o^ + 2 



6 



m 
((-1)^6 



771 



1 



17 \ 1 

-2 log m H h 47 - 

6 J e. 



m 



:d.i2) 
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